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The goal of this note is to state the optimal decay rate for solutions 
of the nonlinear fast diffusion equation and, in self-similar variables, 
the optimal convergence rates to Barenblatt self-similar profiles and 
their generalizations. It relies on the identification of the optimal 
constants in some related Hardy-Poincare inequalities and concludes 
a long series of papers devoted to generalized entropies, functional 
inequalities and rates for nonlinear diffusion equations. 
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Introduction 

The evolution equation 

dr 



[u 



- 1 Vu) = — Au r ' 
m 



[1] 



with m / 1 is a simple example of a nonlinear diffusion equa- 
tion which generalizes the heat equation, and appears in a 
wide number of applications. Solutions differ from the linear 
case in many respects, notably concerning existence, regular- 
ity and large-time behaviour. We consider positive solutions 
u(t, y) of this equation posed for r > and y G R d , d > 1. The 
parameter m can be any real number. The equation makes 
sense even in the limit case m = 0, where u m /m has to be 
replaced by log u, and is formally parabolic for all m G R. No- 
tice that [lj is degenerate at the level u — when m > 1 and 
singular when m < 1. We consider the initial- value problem 
with nonnegative datum u(t = 0,-) = uq G L\ oc (dx), where 
dx denotes Lebesgue's measure on R d . Further assumptions 
on uq are needed and will be specified later. 

The description of the asymptotic behaviour of the solu- 
tions of [JTJl as r — *■ oo is a classical and very active subject. If 
m — 1, the convergence of solutions of the heat equation with 
mo G L\(dx) to the Gaussian kernel (up to a mass factor) is 
a cornerstone of the theory. In the case of equation [TJ with 
m > I, known in the literature as the porous medium equation, 
the study of asymptotic behaviour goes back to [TJ . The result 
extends to the exponents m G (m c , 1) with m c := (d — 2)/d; 
see [2]. In these results the Gaussian kernel is replaced by 
some special self-similar solutions Ud,t known as the Baren- 
blatt solutions (see [3]) given by 



Ud,t(t,u) 



D + 



[2] 



R(r) d V" ' 2d|m— m c \ I R(t) I 

whenever m > m c and ro/1, with 

R(r) := (r + T) , 

where T > and D > are free parameters. To some extent, 
these solutions play the role of the fundamental solution of 
the linear diffusion equations, since lim T ^o Ud,o(t, y) — MS, 
where S is the Dirac delta distribution, and M depends on D. 
Notice that the Barenblatt solutions converge as m — > I to 
the fundamental solution of the heat equation, up to the mass 
factor M. The results of [TJ [2] say that Ud,t also describes 



the large time asymptotics of the solutions of equation [TJ 
as r — > oo provided M = J Rd uo dy is finite, a condition that 
uniquely determines D = D(M). Notice that in the range 
m > m c , solutions of [l] with uq G L\(dx) exist globally in 
time and mass is conserved: J Rd u(t, y) dy — M for any r > 0. 

On the other hand, when m < m c , a natural extension for 
the Barenblatt functions can be achieved by considering the 
same expression ||2|. but a different form for R, that is 



R(t) := (T-r) d (mc ~ m) , 

The parameter T now denotes the extinction time, a new and 
important feature. The limit case m = m c is covered by 



\y\ 2 /d) 



-d/2 



See 



R(t) = e T , U D , T {T,y) = e~ dT (D + e- 
4, 5 for more detailed considerations. 

In this note, we shall focus our attention on the case m < I 
which has been much less studied. In this regime, [TJ is known 
as the fast diffusion equation. We do not even need to assume 
m > 0. We shall summarize and extend a series of recent re- 
sults on the basin of attraction of the family of generalized 
Barenblatt solutions and establish the optimal rates of con- 
vergence of the solutions of [TJ towards a unique attracting 
limit state in that family. To state such a result, it is more con- 
venient to rescale the flow and rewrite [TJ in self-similar vari- 
ables by introducing for m 7^ m c the time-dependent change 
of variables 



t 



loi 



\R(o)J 



and 



1- 



y 



2 d \m — m 



R(t)' ' 

rjd and x 



3] 



with R as above. If m = m c , we take t 
e~ T y/yd. In these new variables, the generalized Barenblatt 
functions UD,T(j,y) are transformed into generalized Baren- 
blatt profiles Vd(x), which are stationary: 



Vd(x 

If u is a solution to [TJ , 

v(t, X 

solves the equation 
dv 



the function 



x G 



[4] 



R(t) u(r,y) 



at 



t>0, x G 



with initial condition v(t = 0,x) — vo(x) := R(0)~ uo(y) 
where x and y are related according to [3] with r = 0. This 
nonlinear Fokker-Planck equation can also be written as 

dv 1 „ m 2 

■ V • (XV) 



at 



1 a n 

— Av 

m 



+ 



1 — 771 



t > 0, x G 



Main results 

Our main result is concerned with the sharp rate at which a 
solution v of the rescaled equation [5j] converges to the gener- 
alized Barenblatt profile Vd given by formula [4[ in the whole 
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range m < 1. Convergence is measured in terms of the relative 
entropy given by the formula 



SM := 



1 



Vh 



D («■ 



V D ) 



dx 



for all m 7^ (modified as mentioned for m = 0). In order 
to get such convergence we need the following assumptions on 
the initial datum vo associated to |f5|: 

(HI) Vd < wo < Vd 1 for some D > Di > 0, 

(H2) if d > 3 and m < m*, (yo — Vd) is integrable for a 

suitable D G [Di,Do]. 

The case m = m* := (d— 4) / (d— 2) will be discussed later. Be- 
sides, if m > m», we define D as the unique value in [D%, Do] 
such that J Rd (vo — Vd) dx — 0. 

Theorem 1. Under the above assumptions, if m < 1 and m 7^ 
m*, t/ie entropy decays according to 



£[v(t,-)] < d 



Vi > 0. 



[61 



27ie sharp decay rate A is equal to the best constant A a>c j > 
in the Hardy-Poincare inequality of Theorem [2] with a := 
l/(m — 1) < 0. Moreover, the constant C > depends only 
on m,d, Do , Di , D and £ [vo] . 

The precise meaning of what sharp rate means will be 
discussed at the end of this paper. As in 4 , we can de- 
duce from Theorem [1] rates of convergence in more standard 
norms, namely, in L q {dx) for q > max{l,d(l — m)/[2(2 — 
m) + d(l — m)]}, or in C k by interpolation. Moreover, by 
undoing the time-dependent change of variables J3J, we can 
also deduce results on the intermediate asymptotics for the 
solution of equation [l] ; to be precise, we can get rates of de- 
cay of u(t, y) — R(r)~ d Ud,t{t, y) as r — > +00 if m G [m c , 1), 
or as r — > T if m G (— 00, m c ). 

It is worth spending some words on the basin of attrac- 
tion of the Barenblatt solutions Ud,t given by J2"]. Such 
profiles have two parameters: D corresponds to the mass 
while T has the meaning of the extinction time of the so- 
lution for m < m c and of a time-delay parameter otherwise. 
Fix T and D, and consider first the case m* < m < 1. The 
basin of attraction of t/o.T contains all solutions correspond- 
ing to data which are trapped between two Barenblatt pro- 
files Ud o ,t(0, •), £^£>i,t(0, ■) for the same value of T and such 
that j Rd [u - Un, T (0,-)]dy = for some D G [D X ,D ]. If 
m < m* the basin of attraction of a Barenblatt solution con- 
tains all solutions corresponding to data which, besides being 
trapped between Ud ,t and Ud 1: t, are integrable perturba- 
tions of Ud,t(0, •). 

Now, let us give an idea of the proof of Theorem [T] 
First assume that D — 1 (this entails no loss of generality). 
On R d , we shall therefore consider the measure dfi a := h a dx, 
where the weight h a is the Barenblatt profile, defined by 
h a (x) := (I + \x\ 2 ) a , with a = l/(m — I) < 0, and study 
on the weighted space L 2 (dfi a ) the operator 

C- a ,d := — hi- a div [ h a V- ] 

which is such that f Rd f (£ a ,d f) dfx a -i = J Rd |V/| 2 d/J, a . This 
operator appears in the linearization of ]5| if, at a formal 
level, we expand 

v(t,x) = h a (x) [l + e f(t,x) /i^" m (x)] 

in terms of e, small, and only keep the first order terms: 



The convergence result of Theorem \T\ follows from the en- 
ergy analysis of this equation based on the Hardy-Poincare 
inequalities that are described below. Let us fix some nota- 
tions. For d > 3, let us define a* := — (d — 2)/2 corresponding 
to m = m»; two other exponents will appear in the analy- 
sis, namely, mi := (d — l)/d with corresponding on = — d, 
and 77i2 := d/(d + 2) with corresponding a?2 = — (d + 2)/2. 
We have m* < m c < 7712 < mi < 1. Similar definitions for 
d = 2 give m„ = —00 so that a„ = 0, as well as m c = 0, 
and mi — 777.2 = 1/2. For the convenience of the reader, a 
table summarizing the key values of the parameter m and the 
corresponding values of a is given in the Appendix. 

Theorem 2. (Sharp Hardy-Poincare inequalities) Let d > 3. For 

any a G (— 00, 0) \ {a*}, there is a positive constant A a ,d such 
that 

Ac, d f |/| 2 dfi a -i < f \Vf\ 2 dfi a MfeH\d^ a ) [7] 

under the additional condition j Rd f d/j, a -i = if a < a*. 
Moreover, the sharp constant A ai d is given by 

( i(d-2 + 2a) 2 if 06 [-^±2,a.)U(a.,0), 

A a , d =i -4a -2d if a G [-d,-^) , 

[ -2a if a G (-00, -d) . 

For d = 2, inequality [7[ /io/ds /or a// a < 0, urai/i t/ie corre- 
sponding values of the best constant A a> 2 = a 2 for a G [—2, 0) 
and A a ,2 ~ —2a for a G (—00, —2). For d = 1, ]7j holds, but 
the values of A a ,i are given by A a ,i = —2a if a < —1/2 and 
A a> i = (a- 1/2) 2 if a G [-1/2,0). 

The Hardy-Poincare inequalities [7[ share many proper- 
ties with Hardy's inequalities, because of homogeneity rea- 
sons. A simple scaling argument indeed shows that 



A Q , d / |/| 2 (D+H 2 ) a - 



x dx< I |V/| 2 (D+ \x[ 



' dx 



holds for any / G H 1 ((D + \x\ 2 ) a dx) and any D > 0, un- 
der the additional conditions J Rd f (D + Ixj 2 )" -1 dx — and 
D > if a < a*. In other words, the optimal constant, A a ,d, 
does not depend on D > and the assumption D = 1 can 
be dropped without consequences. In the limit D — > 0, they 
yield weighted Hardy type inequalities, cf. [BJE]- 

Theorem [2] has been proved in [4] for m < m*. The main 
improvement of this note compared to [5] 2] is that we are 
able to give the value of the sharp constants also in the range 
(777*, 1). These constants are deduced from the spectrum of 
the operator £ a ,d, that we shall study below. 

It is relatively easy to obtain the classical decay rates of 
the linear case in the limit m — * 1 by a careful rescaling such 
that weights become proportional to powers of the modified 
expression (1 + (1 — m) j:rj 2 )~ 1/ ' ( - 1 ~ m ' . In the limit case, we 
obtain the Poincare inequality for the Gaussian weight. As 
for the evolution equation, the time also has to be rescaled by 
a factor (1 — m). We leave the details to the reader. See [9] for 
further considerations on associated functional inequalities. 

A brief historical overview 

The search for sharp decay rates in fast diffusion equations has 
been extremely active over the last three decades. Once plain 
convergence of the suitably rescaled flow towards an asymp- 
totic profile is established (cf. JTJ [2] for m > m c and [U [10] 
for 777 < 777 c ), getting the rates is the next step in the asymp- 
totic analysis. An important progress was achieved by M. 
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Del Pino and J. Dolbeault in pTJ by identifying sharp rates 
of decay for the relative entropy, that had been introduced 
earlier by J. Ralston and W.I. Newman in [121113]. The anal- 
ysis in pj] uses the optimal constants in Gagliardo-Nirenberg 
inequalities, and these constants are computed. J. A. Car- 
rillo and G. Toscani in [14] gave a proof of decay based on 
the entropy/entropy-production method of D. Bakry and M. 
Emery, and established an analogue of the Csiszar-Kullback 
inequality which allows to control the convergence in L 1 (dx), 
in case m > 1. F. Otto then made the link with gradient 
flows with respect to the Wasserstein distance, see |15| . and 
D. Cordero-Erausquin, B. Nazaret and C. Villani gave a proof 
of Gagliardo-Nirenberg inequalities using mass transportation 
techniques in [16] . 

The condition m > (d—l)/d=: mi was definitely a strong 
limitation to these first approaches, except maybe for the 
entropy/entropy-production method. Gagliardo-Nirenberg 
inequalities degenerate into a critical Sobolev inequality for 
m = mi, while the displacement convexity condition requires 
m > mi. It was a puzzling question to understand what 
was going on in the range m c < m < mi, and this has been 
the subject of many contributions. Since one is interested in 
understanding the convergence towards Barenblatt profiles, a 
key issue is the integrability of these profiles and their mo- 
ments, in terms of m. To work with Wasserstein's distance, it 
is crucial to have second moments bounded, which amounts to 
request m > d/(d + 2) =: m 2 for the Barenblatt profiles. The 
contribution of J. Denzler and R. McCann in [171 118] enters 
in this context. Another, weaker, limitation appears when 
one only requires the integrability of the Barenblatt profiles, 
namely m > m c . Notice that the range [m c , 1) is also the 
range for which L l {dx) initial data give rise to solutions which 
preserve the mass and globally exist, see for instance [T]|5]. 

It was therefore natural to investigate the range m € 
(m c , 1) with entropy estimates. This has been done first by 
linearizing around the Barenblatt profiles in .19 , 20 , and then 
a full proof for the nonlinear flow was done by J. A. Carrillo 
and J.L. Vazquez in [21] , A detailed account for these con- 
tributions and their motivations can be found in the survey 
paper [22]. Compared to classical approaches based on com- 
parison, as in the book [5], a major advantage of entropy tech- 
niques is that they combine very well with L^{dx) estimates 
if m > m c , or relative mass estimates otherwise, see [4]. 

The picture for m < m c turns out to be entirely differ- 
ent and more complicated, and it was not considered until 
quite recently. First of all, many classes of solutions vanish 
in finite time, which is a striking property that forces us to 
change the concept of asymptotic behaviour from large-time 
behaviour to behaviour near the extinction time. On the other 
hand, L x (dx) solutions lose mass as time evolves. Moreover, 
the natural extensions of Barenblatt's profiles make sense but 
these profiles have two novel properties: they vanish in finite 
time and they do not have finite mass. 

There is a large variety of possible behaviours and many 
results have been achieved, such as the ones described in [5] 
for data which decay strongly as \x\ — > oo. However, as long 
as one is interested in solutions converging towards Baren- 
blatt profiles in self-similar variables, there were some recent 
results on plain convergence: a paper of P. Daskalopoulos and 
N. Sesum, [10] . using comparison techniques, and in two con- 
tributions involving the authors of this note, using relative 
entropy methods, see [5J[4]. This last approach proceeds fur- 
ther into the description of the convergence by identifying a 
suitable weighted linearization of the relative entropy. In the 
appropriate space, L 2 (d/i Q _i), with the notations of Theo- 
rem [2] it gives rise to an exponential convergence after rescal- 



ing. This justifies the heuristic computation which relates 
Theorems Q] and [2] and allows to identify the sharp rates of 
convergence. The point of this note is to explicitly state and 
prove such rates in the whole range m < 1. 



Relative entropy and linearization 

The strategy developed in [4] is based on the extension of the 
relative entropy of J. Ralston and W.I. Newmann, which can 
be written in terms of w = v/Vn as 



T[w] := 



1 I / r; 
W — 1 I W 

m v 



-1 



1 - m Jmd 

For simplicity, assume m/0. Notice that T\w 
1 



T[w] 



J 

Jjl'I 



V™ dx . 
= S[v]. 

v dx 



Let 



be the generalized relative Fisher information. If v is a solu- 
tion of IJ5J, then 



d — , 



)] = -Z[w(t,-)] \/t>0 



and, as a consequence, lmit^+oo T[w(t, •)] = for all m < 1. 
The method is based on Theorem [5] and uniform estimates 
that relate linear and nonlinear quantities. Following [4] [23] 
we can first estimate from below and above the entropy T in 
terms of its linearization, which appears in [7]] : 



h" 



\f\ 2 Vl~ m dx<2T\w\ < h 2 



f\ 2 Vl~ m dx [9] 



hi(t) ■— inf K dU>(V), h 2 (t) 



where / := (w 

sup R dw(t, ■) and h := max{/i2, 1/hi}. We notice that h(t) — » 1 
as t — > +oo. Similarly, the generalized Fisher information sat- 
isfies the bounds 



\Vf\'V D dx<[l + X(h)]X[w] + Y(h) / |/rV£- 

, 2(2-m) 



dx 
1101 



where /hi < h 



■ =: 1 + X(h) and d(l-m) 
[ {h 2 /hi) 2(2 - m) - 1] < d(l - m) [/i 4(2 ~ m) - 1] =: Y{h). No- 
tice that X(l) — Y(l) — 0. Joining these inequalities with 
the Hardy-Poincare inequality of Theorem [2] gives 

[1 + X(h)] 



2 A Q 



Y(h) 



[11] 



as soon as < h < h* := min{/i > : A aj£ j — Y(h) > 0}. On 
the other hand, uniform relative estimates hold, according to 
[23] , formula (5.33): for some C = C(d, m, D, Do, D i), 

< h - 1 < QFd+2-(r+i) m _ [ 12 ] 

Summarizing, we end up with a system of nonlinear differen- 
tial inequalities, with h as above and, at least for any t > t*, 
t» > large enough, 

A a , d - Y(h) 



d ^.r , 



< 



*>(*,-)]• 



13] 



[1 + X(h)] h 2 
Gronwall type estimates then show that 

limsup e 2 A "' d t J r [w{t, •)] < +oo . 

t — >oo 

This completes the proof of Theorem [T] for The adap- 

tation to the logarithmic nonlinearity is left to the reader. 
Results in [J are improved in two ways: a time-dependent es- 
timate of h is used in place of h(0), and the precise expression 
of the rate is established. One can actually get a slightly more 
precise estimate by coupling [1121] and [1131] . 
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Corollary 3. Under the assumptions of Theorem [TJ if h(0) < 
h,, then T[w(t,-)] < G(t, h(0), F[w(0, ■)]) for any t > 0, 
where G is the unique solution of the nonlinear ODE 



dG 



= -2; 



Y(h) 



G with h = 1 + CG d + 2 -<- d + 1 ^ 



dt [1 + X{h)\ h 2 

and initial condition G(0) = J-[w(0, ■)]. 

Operator equivalence. The spectrum of C* a ,d 
An important point of this note is the computation of the 
spectrum of C a< d for any a < 0. This spectrum was only par- 
tially understood in 4, 8]. In particular, the existence of a 
spectral gap was established for all a 7^ a* = (2 — d)/2, but 
its value was not stated for all values of a. 

J. Denzler and R.J. McCann in [17LI18] formally linearized 
the fast diffusion flow (considered as a gradient flow of the en- 
tropy with respect to the Wasserstein distance) in the frame- 
work of mass transportation, in order to guess the asymptotic 
behaviour of the solutions of flTJ. This leads to a different 
functional setting, with a different linearized operator, TL a ,d- 
They performed the detailed analysis of its spectrum for all 
m G (m c , 1), but the justification of the nonlinear asymptotics 
could not be completed due to the difficulties of the functional 
setting, especially in the very fast diffusion range. 

Our approach is based on relative entropy estimates and 
the Hardy- Poincare inequalities of Theorem [2] The asymp- 
totics have readily been justified in [4|. The operator C a ,d = 



-h\~ a div [h a V- ] can be initially defined on 



). To con- 



struct a self-adjoint extension of such an operator, one can 
consider the quadratic form / 1— > (f,C a ,df), where (•, •) 
denotes the scalar product on L 2 (dfi a -i). Standard results 
show that such a quadratic form is closable, so that its clo- 
sure defines a unique self-adjoint operator, its Friedrich's ex- 
tension, still denoted by the same symbol for brevity. The 
operator H ay d is different: it is obtained by taking the opera- 
tor closure of C a ,d, initially defined on D(R d ), in the Hilbert 
space H 1 '*(dfi a ) := {/ G L 2 (d/^ a _i) : J Rd |V/| 2 dfi a < 00 and 
J R d = if a < a,}, so that 

H a , d f ■= hi- a V ■ [h a V(hi- a V ■ (h a V/))] . 

Because of the Hardy-Poincare inequality, J Rd |V/| 2 djx a de- 
fines a norm. Denote by ( ■ , ■ } the corresponding scalar prod- 
uct and notice that ( / , g ) = ( / , C a ,d g ) for any /, g G V(R d ) . 

Proposition 4. The operator C a ,d on L 2 (dfi a ~i) has the same 
spectrum as the operators Ti. a ,d on H 1 '*(dfi a ). 

The proof is based on the construction of a suitable 
unitary operator U : H 1 '*(djj, a ) — + L 2 (djj, a -i), such that 
U TLa,dU~ l = C a ,d- We claim that U = C a ,d is the re- 
quested unitary operator. By definition, D(R d ) is a form 
core of C a ,d, and as a consequence, the identity has to be 
established only for functions / G £>(R d ). Since [\U f[\ 2 = 
(/, C a ,df) = J Rd |V/| 2 d/i a , we get 

{UHa^U- 1 / ,g) = (Ha^U- 1 /, U-'g) 

= [C^dU^f ^C^dU^g) = (Uf,Ug) = (f,g) 

where we have used the properties U* = U -1 and U 2 = C a ,d- 
This unitary equivalence between £ QjC j and 7i a ,d implies the 
identity of their spectra. 

We may now proceed with the presentation of the actual 
values of the spectrum by extending the results of [18] . Ac- 
cording to [24] . the spectrum of the Laplace-Beltrami operator 



on S d 1 is described by 

-A sd -iYi fl = l(e + d-2)Y ifl 

with^O, I, 2, . . . and fi — 1, 2, . . . M t := (d+t ~*f d ( _ d +f l - 2) 
with the convention Mo = I, and Afi = 1 if d = 1. Using 
spherical coordinates and separation of variables, the discrete 
spectrum of C a ,d is therefore made of the values of A for which 



" 1 
v + 



d—l 1 2 a r 
r ' 1 + r 2 



V + 



X 

T+r 2 



i (l+d-2) 

— 71 — 







[141 



has a solution on 

of variables v(r) = r l w(—r 2 ) allows to express w in terms of 
the hypergeometric function 2-F\(a, b, c; 2) with c = t + d/2, 
a + b+1 = £ + a + d/2 and a b = (2 £ a + A)/4, as the solution 
for s = — r 2 of 

s (1 — s) y" + [c — (a + b + I) s] y' — aby = , 

see [25]. Based on [HJ 18 , we can state the following result. 

Proposition 5. The bottom of the continuous spectrum of the 
operator C a ^d on L 2 (dfi a -i) is A^f := l(d + 2a-2) 2 . More- 
over, Ca.d has some discrete spectrum only for m > m2 = 
d/(d + 2). For d > 2, the discrete spectrum is made of the 
eigenvalues 



A f 



-2a(£ + 2fc)-4fc [k + £ + 



d 



[151 



with I, k = 0, 1, . . . provided (£, k) (0, 0) and £ + 2k - 1 < 
-(rf + 2a)/2. If d — 1, the discrete spectrum is made of the 
eigenvalues Afe = k (1 — 2 a — k) with k G N PI [1, 1/2 — a] . 

Using Persson's characterization of the continuous spec- 
trum, see |26ll8j. one can indeed prove that A^™' i s the optimal 
constant in the following inequality: for any / G T>(M> d \ {0}), 



A c „ 



l/i 



.|2(a-l) 



dx < 



|V/| 2 |a:| 2Q d:r. 




■ d 

Fig. 1. Spectrum of C a> d as 3 function of a, for d = 5. 
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The condition that the solution of | I14| | is in the domain of 
Ca.d determines the eigenvalues. A more complete discussion 
of this topic can be found in |18| . which justifies the expression 
of the discrete spectrum. 

Since a = l/(m — 1), we may notice that for d > 2, a = — d 
(corresponding to —2a = Aio = Aoi = —4a — 2d) and a = 
-(d+2)/2 (corresponding to A i = A[j ont := | (d+2a-2) 2 ) re- 
spectively mean m = mi = (d— l)/d and m = mi = d/(d+2). 

The above spectral results hold exactly in the same form 
when d = 2, see |18| . Notice in particular that A„°3L2 — q2 
so that there is no equivalent of m, for d = 2. With the no- 
tations of Theorem [2] a* — 0. All results of Theorem [1] hold 
true under the sole assumption (HI). 

In dimension d = 1, the spectral results are different, see 
[18] . The discrete spectrum is nonempty whenever a < —1/2, 
that is m > — 1. 



The critical case 

Since the spectral gap of C a ^d tends to zero as m — t m* , 
the previous strategy fails when m = m* and one might ex- 
pect a slower decay to equilibrium, sometimes referred as slow 
asymptotics. The following result has been proved in [23] . 

Theorem 6. Assume that d > 3, let v be a solution of ]5j with 
m = m„, and suppose that (H1)-(H2) hold. If \vo — Vd\ is 
bounded a.e. by a radial L (dx) function, then there exists a 
positive constant C* such that 

£[v{t,-)]<CT 1/2 Vt>0, [16] 

where C* depends only on m, d, Do, D%, D and £\vq\. 

Rates of convergence in L q (dx), q G (l,oo] follow. Notice 
that in dimension d = 3 and 4, we have respectively m, = — 1 
and m, = 0. In the last case, Theorem [6] applies to the loga- 
rithmic diffusion. 

The proof relies on identifying first the asymptotics of the 
linearized evolution. In this case, the bottom of the contin- 
uous spectrum of C a ,,d is zero. This difficulty is overcome 
by noticing that the operator C a ,.d on R d can be identified 
with the Laplace-Beltrami operator for a suitable conformally 
flat metric on M d , having positive Ricci curvature. Then the 
on-diagonal heat kernel of the linearized generator behaves 
like t~ d/2 for small t and like t~ 1/2 for large t. The Hardy- 
Poincare inequality is replaced by a weighted Nash inequality: 
there exists a positive continuous and monotone function A/" 
on R + such that for any nonnegative smooth function / with 
M = J Rd fdfi-d/2 (recall that a, — 1 = — d/2), 




The function M behaves as follows: lim s _ + s 1//3 A/"(s) = 
ci > and lim^oo s" d/(d+2) A'(s) = c 2 > 0. Only the first 
limit matters for the asymptotic behaviour. Up to technicali- 
ties, inequality | |13| | is replaced by (T[w(t, -)]) 3 < KX[w(t, •)] 
for some K > 0, t > to large enough, which allows to complete 
the proof. 

Faster convergence 

A very natural issue is the question of improving the rates of 
convergence by imposing restrictions on the initial data. Re- 
sults of this nature have been observed in [TS] in case of radi- 
ally symmetric solutions, and are carefully commented in |18j . 
By locating the center of mass at zero, we are able to give an 
answer, which amounts to kill the Aio mode, whose eigenspace 



is generated by x i— > Xi, i = 1, 2. . . d. This is an improvement 
compared to the first result in this direction, which has been 
obtained by R. McCann and D. Slepcev in [27], since we obtain 
an improved sharp rate of convergence of the solution of [5[ , 
as a consequence of the following improved Hardy-Poincare 
inequality. 

Lemma 1. Let A QjC j := —4a — 2d if a < —d and A. a ,d '■= A^°Jf 
if a G [— d, — d/2). If d > 2, for any a G (— oo,— d), we have 

A a , d f |/i 2 d Ma _i< f |V/| 2 d/x Q V/6fl'(W 

under the conditions f Rd /d/i a _i = and J Rd x f d/^ Q _i = 0. 

The constant A a> d is sharp. 

This covers the range m G (mi, 1) with mi = (d — l)/d. 

Theorem 7. Assume that m G (mi, 1), d > 3. Under As- 
sumption (HI), if v is a solution of [5|] with initial da- 
tum vo such that J Rd x vo dx = and if D is chosen so that 

Jfid( v o — Vd) dx = 0, then there exists a positive constant C 
depending only on m, d, Do, D\, D and £[vo] such that the rel- 
ative entropy decays like 

e[v(t,-)]<Ce- Aa - dt Vt>0. 
A variational approach of sharpness 

Recall that (d — 2)/d = m c < mi = (d — l)/d. The entropy 
/ entropy production inequality obtained in [11] in the range 
m G [mi, 1) can be written as T < ~I and it is known to be 
sharp as a consequence of the optimality case in Gagliardo- 
Nirenberg inequalities. Moreover, equality is achieved if and 
only if v = Vd- The inequality has been extended in [21] to 
the range m G (m c , 1) using the Bakry-Emery method, with 
the same constant 1/2, and again equality is achieved if and 
only if v = Vd, but sharpness of 1/2 is not as straightforward 
for m G (m c ,mi) as it is for m G [mi, 1). The question of 
the optimality of the constant can be reformulated as a vari- 
ational problem, namely to identify the value of the positive 
constant 



where the infimum is taken over the set of all functions such 
that v G 0(R d ) and J Rd v dx = M. Rephrasing the sharpness 
results, we know that C — 2 if m G (mi, 1) and C > 2 if 
m G (m c , mi). By taking v n = Vd (1 + ^ / V D ~ m ) and letting 
n — > oo, we get 

T[v n ] _ J Rd \Vf\ 2 Vpdx 
n^£[v n ] f Rd \f\*V 2 - m dx- 

With the optimal choice for /, the above limit is less or equal 
than 2. Since we already know that C > 2, this shows that 
C — 2 for any m > m c . It is quite enlightening to observe 
that optimality in the quotient gives rise to indetermination 
since both numerator and denominator are equal to zero when 
v = Vd- This also explains why it is the first order correction 
which determines the value of C, and, as a consequence, why 
the optimal constant, C = 2, is determined by the linearized 
problem. 

When m < m c , the variational approach is less clear since 
the problem has to be constrained by a uniform estimate. 
Proving that any minimizing sequence (v n ) n evi is such that 
Vti/Vd — 1 converges, up to a rescaling factor, to a function / 
associated to the Hardy-Poincare inequalities would be a sig- 
nificant step, except that one has to deal with compactness 
issues, test functions associated to the continuous spectrum 
and a uniform constraint. 
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Sharp rates of convergence and conjectures 

In Theorem[T] we have obtained that the rate exp(— A QiC j t) is 
sharp. The precise meaning of this claim is that 

v • * - * T M 

A a ,d = limint mf , 

h->0 + J- [W\ 

where the infimum is taken on the set Sh of smooth, nonneg- 
ative bounded functions w such that \\w — l||£°°(<fa) < h and 
such that J R d(w — 1) Vd da; is zero if d = 1,2 and m < 1, 
or if d > 3 and m» < m < 1, and it is finite if d > 3 and 
m < m*. Since, for a solution x) = w(t,x)VD(x) of 
|j5j, [8j] holds, by sharp rate we mean the best possible rate, 
which is uniform in t > 0. In other words, for any A > A a> d, 
one can find some initial datum in Sh such that the estimate 
J-[w(t, ■)] < JF[to(0, ■)] exp(— \t) is wrong for some t > 0. We 
did not prove that the rate exp(— A a ,d t) is globally sharp in the 
sense that for some special initial data, !F[w(t, •)] decays ex- 
actly at this rate, or that lim inft^oo exp(A a ,d t) J-\w(t, ■)] > 0, 
which is possibly less restrictive. 

However, if m £ (mi,l), mi = (d — l)/d, then 
exp(— A a .d t) is also a globally sharp rate, in the sense that 
the solution with initial datum uo(x) = Vd(x + xo) for any 
xo £ K \ {0} is such that !F[w(t, ■)] decays exactly like 
exp(— A. a ,d t). This formally answers the dilation-persistence 
conjecture as formulated in [18]. The question is still open 
when m < mi. 

Another interesting issue is to understand if improved 
rates, that is rates of the order of exp(— A« t) with (£, k) ^ 



(0, 0), (0, 1), (1, 0) are sharp or globally sharp under additional 
moment-like conditions on the initial data. It is also open to 
decide whether exp(— K a ^t) is sharp or globally sharp under 
the extra condition J Rd xvodx = 0. 



Appendix. A table of correspondence 

For the convenience of the reader, a table of definitions of the 
key values of m when d > 3 is provided below with the cor- 
respondence for the values of a — l/(m — 1). This note is 
restricted to the case m G (— oo, 1), that is a £ {—oo, 0). 



m = 


— oo 


m„ 


m c 


m 2 


mi 


1 


m 


— oo 


d-A 

d-2 


d-2 
d 


d 

d+2 


d-1 
d 


1 


■ = 





d-2 
2 


d 

2 


d+'l 

2 


-d 


— oo 
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